Identification of Nuclear Relaxation Processes in a Gapped Quantum Magnet: 
X H NMR in the S = \ Heisenberg Ladder Cu 2 (C 5 H 12 N 2 ) 2 Cl4 
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The 1 H hyperfine shift K and NMR relaxation rate Tj" 1 have been measured as a function of 
temperature in the S = 1/2 Heisenberg antiferromagnetic ladder Cu2(C5Hi2N2)2Ci4. The presence 
of a spin gap A ~ J± — Jm in this strongly coupled ladder (Jm < J±) is supported by the K and 
Tj" 1 results. By comparing T^ 1 at two different 1 H sites, we infer the evolution of the spectral 
functions S z {q,ui n ) and S±(q,u> n )- When the gap is significantly reduced by the magnetic field, two 
different channels of nuclear relaxation, specific to gapped antiferromagnets, are identified and are 
in agreement with theoretical predictions. 



Several classes of one-dimensional Heisenberg antifer- 
romagnets (HAF) are known to exhibit a spin-gap at low 
temperature. For example, integer-spin chains Q have 
a non-magnetic " spin- liquid" ground state (singlet) sep- 
arated from a branch of triplet excitations by an energy 
gap A. A spin- liquid ground state also exists in spin- 
ladders, built by coupling an even number of S = 1/2 
HAF chains with an antiferromagnetic transverse ex- 
change Jl [@,|] . At low energies, many physical prop- 
erties are dominated by the singlet-triplet gap and do 
not depend on the underlying dynamical quantum pro- 
cesses stabilizing the ground state. For example, thermo- 
dynamic quantities (susceptibility, specific heat) are very 
similar in a number of gapped one-dimensional HAF. 

In this letter, the low-energy dynamical processes dom- 
inating the : H spin-lattice relaxation rate (1/Ti) of an 
organic spin-ladder (Cu2(C5Hi2N2)2Ci4) are unambigu- 
ously identified by comparing the T\ measurements at 
two proton sites. They coincide precisely with the pro- 
cesses proposed by Sagi and Affleck [[| for Haldane sys- 
tems. This experimental evidence supports the idea pro- 
posed by Sachdev and coworkers j^] that spectral func- 
tions S z< ±(q,w) are, at low energies (u -C A), common 
to all gapped one-dimensional HAF. 

In Cu2(C 5 Hi2N 2 )2Cl4 [§, the Cu 2 + (S=l/2) ions are 
coupled antiferromagnetically in well isolated ladders |t]] 
(see Fig. |l|). The exchange parameters along the rungs 
( J_l) and the legs (Jm ) of the ladder are isotropic and have 
been accurately measured to be Jj_ = 13.2 K and Jn = 
2.4 K ^-jT^]. In many respects, this material is a model 
system in which theoretical predictions for Heisenberg 
ladders in the strong coupling limit (J±/J\\ = 5.5 ^> 1) 
can be tested. 

The 1 H NMR measurements were carried out by pulsed 
spin-echo techniques on five single crystals (typically 
1 x 1 x 0.05 mm 3 ), oriented with their b axis (perpendic- 



ular to the chains axes) along the applied field Hq = 5.6 
Tesla. A typical spectrum for the proton resonance shows 
a number of partially resolved lines (Fig. 2a), indicating 
a variety of local fields among the 24 inequivalent *H 
sites. In the following, we focus on the lines labelled (I) 
and (II), as their extreme position in the spectrum al- 
lows their study on a wide temperature (T) range flTT| |. 
With the field along 6, the magnetic hyperfine shift 
of the proton resonance is related to the uniform spin 
susceptibility xo = X*(<7 = 0' w = 0) a ^ the nuclear site i, 
by 
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where A^b is the hyperfine coupling constant and a the 
chemical shift. The shift of the two lines, plotted in 
Fig. 2b, are opposite in sign but follows the same In- 
dependence as xoj that is, a high temperature Curie- 
Weiss behavior followed by an exponential drop below a 
rounded maximum at T x m^ ~ 8 K, in complete agree- 
ment with previous susceptibility measurements p|,^0[. 
Since K is proportional to the susceptibility \o measured 
at 5 T (Inset to Fig. 2b), the hyperfine couplings on 
both sites can be estimated: a\^ = 



and A ( b [ I] = -2.6 ± 0.50 kOe |1 



2.95 ± 0.40 kOe 
The largest contri- 
bution to A comes from the dipolar field on the j-th 
nucleus created by the surrounding electronic spins, i.e. 

Aj cx -|7 e |7„?i 2 J2i i 1 ~ 3cos2 %)/ r rr wnere % is tne 
angle between and Hq. Given the atomic positions, 
it is straightforward to compute the dipolar field at each 
4 H site (a reliable result is obtained by summing over 
5-6 neighboring Cu spins). The total NMR spectrum is 
well-reproduced in this way. It is therefore possible to 
assign the NMR lines to specific proton sites: line (II) 
is ascribed to protons H2 involved in the superexchangc 
Jm (See Fig. |l|). The line (I) is attributed to protons 
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H20 and H23 at the outer edges of the ethyl groups. It 
must be stressed that the uncertainty in the site labelling 
could only result in adding the protons H14 and H4 to 
the lines (I) and (II), respectively. This essentially does 
not affect our analysis of the nuclear relaxation. 

In a magnetic field, the triplet excitations split into 
three branches. In the strong coupling limit (J\\/Jj_ <C 
1), the lowest branch is separated from the singlet ground 
state by an effective gap Ah, = A — h, where h — g^sHo is 
the Zeeman energy. When the temperature is small com- 
pared to A/j, interactions between excitations are negli- 
gible and the lowest branch dominates the temperature 
dependence of the susceptibility Jl3| 
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A low temperature fit of to Eq. ([[]) and (^) gives an 
effective gap of Ah ~ 3 K in 5.6 Tesla. This is very close 
to the value expected taking the zero-field gap A = 10.8 
K inferred from susceptibility and high field magneti- 
zation measurement Q reduced by the Zeeman energy 
h = 7.6 K. Thus, the measurements of K and Xo are 
fully consistent with a spin gap A^ = ( J± — Jm ) — h 
between the singlet and triplet states of the Heisenberg 
ladder with strong rungs. 

We now discuss the dynamical properties, as probed 
by the nuclear spin-lattice relaxation rate 1 fT\ . The re- 
covery of the nuclear magnetization is always a single 
exponential at all temperatures. As shown in Fig. 3a, 
both sites display qualitatively the same T-dependence, 
that is, 1/Ti tends to be constant in the paramagnetic 
region, and crosses over to an activated behavior at low 
T . There are, however, striking differences between the 
two lines: (a) at high-T the values of T\ differ by one 
order of magnitude, (b) at low-T the gap values differ 
by a factor of 2. Indeed, assuming an activated behavior 
1/Ti oc exp(-A e// /fc B T), A eff ~ 3.4 ± 0.2 K for the 
line (I), close to the value Ah = 3.0 K deduced from the 
shift, but A eff ~ 6.8 ± 0.2 K for the line (II). 

In spin systems, the temporal fluctuations of the elec- 
tronic spins make the nuclear polarization relax in a time 
Ti related to the spectral densities S z ±(q,uj) of the two- 
spin correlation functions, through ]l4| 



1 



With UJn 
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[F z S z (q,LU n )+F ± S ± (q,uj n )}, (3) 



the nuclear Larmor frequency, and 



S z ,x(q,u n )= dte^dS^+iq^S^i-q^)}). (4) 



In general, any quantitative analysis of the relaxation 
requires the knowledge of the hyperfine "form" factors 
F(q) [jl5| in addition to a model for S Z} ±(q, uo). We first 
discuss the structure factors S z ,±(q,cu). 



Single magnon processes, which require an energy 
greater or equal to the gap, cannot contribute to the nu- 
clear relaxation which involves negligible energy transfers 
huj n ~ mK. Two- or three-magnon scattering processes 
are then required |ll| . More specifically, Sagi and Affleck 
[pj have recently analyzed the possible nuclear relaxation 
processes for Haldane chains in magnetic fields. Since the 
low energy excitations of S = 1/2 ladders and integer- 
spin chains are qualitatively similar Jl7| , it is natural to 
consider the same processes here. Following their argu- 
ments, the nuclear spins can exchange energy through 
three different channels [Q: 

- (i) "intrabranch" transitions involve two magnons 
within the same branch (i.e. with the same S z eigen- 
value). At low T, these processes have a maximum 
probability near the minimum at k = n of the low- 
est branch of the triplet (Fig. ||), implying a momen- 
tum transfer Afc = q ~ (forward scattering). For 
T <C h, the g-integrated spectral density is expected 
to follow the thermal occupation of the lowest energy 
triplets Sl ntra (uj n ) cx exp (-(A - h)/k B T). 

- (ii) "interbranch" transitions (Fig. |]), i.e. transitions 
from a state in a magnon branch m to a state with mil 
(S± operators). Since the Zeeman splitting at 5.6 Tesla is 
larger than the magnon bandwidth (~ 5.5 K @), these 
processes can only occur because of the finite damping 
of each level and are expected to be weak. Furthermore, 
only large momentum transfers q ~ tt (backward scatter- 
ing) remain at large Zeeman splitting. One infers from 
§ that Sf ter (u n ) oc exp(-A//c B T). 

- (Hi) "staggered" processes: when Hq approaches the 
critical field h c \ « A, one-magnon excitations (S± op- 
erators, q ~ 7r) become increasingly relevant. At fi- 
nite T in the gapped phase, interactions between ex- 
citations, or equivalently finite damping, generate non- 
vanishing matrix elements for such transitions: this re- 
laxation mechanism involves three-magnon (or higher or- 
der) processes, and its temperature dependence follows 
the square of the thermal population in the lowest triplet 
state, Sf a "{uJ n ) oc exp(-2(A - h)/{k B T)). 

The above discussion shows that: (1) the Boltzmann 
factor is more favorable to intrabranch processes (Ah ~ 3 
K); these will dominate the staggered transitions (2A^ ~ 
6 K), while interbranch ones, if any, are essentially neg- 
ligible (A ~ 10 K). (2) The low-T nuclear relaxation is 
only driven by two terms: S z (q ~ 0,w„) for intrabranch 
transitions and S±(q ~ it,Lu n ) for staggered transitions. 
Accordingly, we write Eq. ([|) in a simplified form: 

l/T 1 ocF z (0)Sz(q = 0,u n ) + F x (7r)Sx(q = TT,u n ). (5) 

Hence, the two behaviors 1/Ti oc exp(— Ah/k^T) 
for the line (I) and \jT x oc exp(-2A^/fc B T) for 
(II) can only come from the temperature dependence 
of S z (q = 0,w„) oc exp(-A h /k B T) while S±(q = 
ir,uj n ) oc exp(— 2Ah/k B T). In other words, the ratio of 
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F z (0)/F±(n) for lines (I) and (II) are such that, at low 
temperatures, only one of the exponential terms domi- 
nates the relaxation: obviously, S z (q = 0, w„) component 
is favored for the line (I), and S±(q — n,uj n ) one for the 
line (II). 

This result is, to our knowledge, the first experimen- 
tal identification of specific nuclear relaxation channels in 
a gapped antiferromagnet, a result in remarkable agree- 
ment with the work of Sagi and Affleck. Another support 
to this theory is that the lowest gap Ah — 3 K is also 
the value seen in the susceptibility. Furthermore, the 
observation of the staggered contribution in the gapped 
phase proves that interactions between fermionic-like ex- 
citations are significant in this system. This conclusion 
was already drawn from magnetization measurements || . 

A nice feature of this study is that the T\ data for 
the two lines provide a set of two independent equations 
(i.e. Eq. || for each line). Since F z (q) and F±(q) can 
be computed for each 1 H site, the spectral functions S z 
and S± can, in principle, be extracted separately. Here, 
we found that F±(q = it) is indeed five times larger than 
F z (q — 0) for the line (II), while both terms are com- 
parable for the line (I). S± oc exp(— 2Ah/k-QT) is thus 
overweighted for the line (II) explaining why T[ J decays 
with an activation energy 2 Ah- On the other hand, T/ 
is predominantly sensitive to the smallest gap, A^, gen- 
erated by S z (q — 0,uj n ). However, one must realize that 
the calculation of the form factors is subject to several 
uncertainties: any error in atomic positions is amplified 
(F(rij) oc r~j 6 ), the spatial extension of Cu 2+ orbitals 
may play an important role |l9| for the protons H2 (line 
II) which are in the superexchange pathway and closer 
to the Cu ion than those of line (I). Indeed, the ex- 
tracted values of S z are slightly negative suggesting that 
the value of Fj 1 has been underestimated in the calcu- 
lation. In fact, the pure Heisenberg paramagnetic limit 
Sz(<l,w n ) — ^S±(q,uj n ) should be recovered when T is 
large compared to J± t \\ and h. F± can be rescaled to 
a value satisfying this limit at T = 30 K, where the ob- 
served value of T\ for the line (I) is within 10 % of the 
paramagnetic limit calculated by Moriya Q]. In any 
event, this rescaling does not affect the gap parameters 
extracted from the low-T behaviour. As shown in Fig. 
3c, S± experiences a gap A ~ 6.8 ± 0.2 K twice as large 
as in S z (3.4 ±0.2 K). 

This analysis of the nuclear relax- 
ation in Cu2(C5Hi2N 2 )2Cl4, a strongly coupled ladder, 
settles an ongoing controversy: the gap values derived 
from static (xo) an d dynamic S z (q ~ 0, u> n ) measure- 
ments are here identical, in complete agreement with the 
predictions of Troyer et al. |b| and Sagi and Affleck [Q . 
This strongly contrasts with the experimental observa- 
tions in inorganic ladders and in some Haldane 
chains ^,^3|. In these materials, the different tempera- 
ture dependence observed in the dynamics may be due to 
a second minimum in the dispersion relation. For exam- 



ple, it should be the case in SrCu203 if J± < J|| p4] , p5| . 
Low-lying excitations near k = such that Ak=o ~ Ak=n 
may open up relaxation channels involving large-g inter- 
branch transitions and would lead to a higher effective 
gap in T\ measurements. Other explanations have been 
proposed in the limits T <C A h @ and J± < Jy 

In conclusion, 1 H NMR experiments demonstrate 
that the effective gap of the S = 1/2 HAF ladder 
Cu2(C5Hi 2 N 2 )2Cl4 in a magnetic field h is A^ ~ J± — 
J|| — h. The nuclear relaxation can be quantitatively un- 
derstood in the framework of the theory of Sagi and Af- 
fleck, retaining only " intrabranch" and "staggered" pro- 
cesses. More generally, the processes identified in this 
work should be generic to many gapped HAF chains in 
a magnetic field. 
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tionne aux Universites J. Fourier et INPG Grenoble I. 
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FIG. 1. Schematic structure of Cu2(CsHi2N2)2Cl4, with 
the exchange parameters determined in Ref. ||. The labelled 
protons contribute to the two NMR lines used in this work to 
probe different dynamical functions (S z and S±, see text). 

FIG. 2. (a): 1 H NMR spectra at fixed frequency 
fo = 239.112 MHz. The arrows indicate the two lines studied 
(their different amplitudes are due to different excitation con- 
ditions). The large central peak comes from protons in the 
NMR probe, (b): magnetic hyperfine shift for the lines (7) 
and (II) and susceptibility at 5 Tesla; the dashed line is a fit 
where the hyperfine coupling is the only adjustable parameter 
and temperature dependence of the susceptibility is given by 
Eq. 5 of Ref. §, with J± = 13.2 K and Jy = 2.5 K. Inset: 
shift data vs. susceptibility, with T as an implicit parameter. 

FIG. 3. X H spin-lattice relaxation rate 1 /Ti for lines (I) and 
(II) as a function of T (a) and 1/T (b). (c): S z (q = 0,u} n ) and 
S±(q — n,u n ) correlation functions derived from Ti results 
(see text for details). 

FIG. 4. Schematic picture of the two-magnon scattering 
processes relevant to the nuclear relaxation in a system with 
singlet to triplet gap, in a magnetic field Ho Q|. In this ex- 
periment, the Zeeman splitting g^sH ~ 7.6 K is larger than 
the magnon bandwidth (~5.5 K |18[|). 
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